Detection of the Meissner Effect with a Diamond Magnetometer 
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We examine the possibility of probing superconductivity effects in metal nanoclusters via dia- 
mond magnetometry. Metal nanoclusters have been proposed as constitutive elements of high T c 
superconducting nanostructured materials. Magnetometry based on the detection of spin-selective 
fluorescence of nitrogen- vacancy (NV) centers in diamond is capable of nanoscale spatial resolution 
and could be used as a tool for investigating the properties of single or multiple clusters interacting 
among each other or with a surface. We have carried out sensitivity estimates and experiments to 
understand how these magnetometers could be used in such a situation. We studied the behavior of 
the sensor as a function of temperature and detected the flux exclusion effect in a superconductor by 
monitoring the [111] -orient at ion fine structure spectrum in high-NV density diamond. Our results 
show that phase transitions can be ascertained in a bulk superconductor with this technique while 
the principal zero field splitting parameter (D) demonstrates a temperature dependence that may 
require compensation schemes. 

PACS numbers: 36.40.-c, 74.25.Nf, 74.25.Ha, 74.25.Dw, 74.25.-q, 75.20.-g 



Superconductivity occurs in a wide variety of mate- 
rials when the temperature drops below a critical tem- 
perature T c and is accompanied by diamagnetism and 
a drop in resistance to electric current. The super- 
conducting transitions observed to date are all below 
135 K under normal atmospheric pressure and they are 
mostly associated with bulk systems. In recent years, 
the search for room-temperature superconductivity has 
been approached from several different angles [l| , includ- 
ing searching for systems with negative dielectric func- 
tions Q , synthesizing crystals with specific phonon spec- 
tra [H, creating materials with a pore structure [4], and 
studying the superconducting state at interfaces in epi- 
taxial heterostructures [5[. 

Metal nanoclusters, which can exhibit superconduc- 
tivity effects 0-0? a ^ so f° rm a novel family and could 
potentially develop into room-temperature superconduc- 
tors. Size-selected metal nanoclusters, in the case of cer- 
tain metals, are known to exhibit electronic shell struc- 
ture [6-8]. Kresin and Ovchinnikov recently proposed 
that a superconducting transition can occur for some 
metal nanoclusters with 10 2 -10 3 valence electrons [ol. [Toj . 
The transitions are expected to occur near electronic shell 
closings. Under special circumstances, the shell structure 
is expected to strengthen pairing correlations, and leads 
to elevated values for the critical temperature for a transi- 
tion to the superconducting state. The superconducting 
state of small nanoclusters is directly related to the phe- 
nomenon of pair correlation. Recently, Cao et al. have 
presented experimental measurements of a jump in the 
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heat capacity of Al^ and AlJ 7 clusters, which is strong 
evidence of a superconducting transition in isolated clus- 
ters [ll|. The transition temperature was T c w 200 K. A 
number of interesting predictions for other clusters have 
been made in 0, [To| . 

The clusters' exceedingly small sizes (1 nm or less) 
make them difficult to study in the laboratory, as they 
are inaccessible to a large number of common spectro- 
scopies and methods for detecting superconductivity. It 
is important to have additional means to probe such clus- 
ters at the individual level and in populations, so that 
cluster-cluster or cluster-surface interactions can be stud- 
ied. Recent developments in solid-state magnetometers 
such as anisotropic magnetoresistive sensors [H, [l3| or 
those based on diamonds with nitrogen-vacancy (NV) 
color centers [l34l6j have opened up new possibilities for 
sensitive magnetic measurements of material properties 
over a broad range of temperatures with high spatial res- 
olution. "NV-diamond" magnetometers have been devel- 
oped which can perform measurements at the nanoscale, 
potentially, with single-spin detection [14], ll6| . 

In this article, we look at the feasibility of performing 
magnetic measurements near the surface of a metal clus- 
ter and bulk superconductor to detect the phase transi- 
tion through the Meissner effect or flux exclusion. There 
are three common approaches to finding the critical tem- 
perature in superconductors. These consist of measure- 
ments of magnetic susceptibility but also electrical re- 
sistivity and specific- heat. Resistivity measurements at 
the nanoscale are challenging. Very often, more than one 
type of experimental measurement is needed to convinc- 
ingly prove the existence and probe the nature of the 
superconducting state. For heat capacity measurements 
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in individual isolated clusters we note the recent work 
of Jarrold et al. [I?], G3 • Susceptibility measurements re- 
port changes in the magnetic moment below T c by way of 
a susceptometer instrument. The measurement directly 
probes the Meissner effect in which a transition to dia- 
magnetism occurs. 

Diamagnetism in metal clusters is different from bulk 
metals. Clusters containing a magic number of electrons 
(i.e. filled shell) are diamagnetic regardless of the pairing 
correlation because of the spherical symmetry of the elec- 
tronic shell structure [9]. This is an important deviation 
from bulk superconductors which indicates there would 
be no detectable Meissner effect. For clusters with in- 
complete shells, however, there is orbital paramagnetism 
at high temperatures [9j and the transition to the su- 
perconducting state at T=T C should be accompanied by 
the paramagnetic-diamagnetic transition in the cluster 
moment. This suggests there may be a detectable Meiss- 
ner effect. On the other hand, the London penetration 
depth for bulk superconductors (typical values are be- 
tween 50 nm and 500 nm), is generally much larger than 
the cluster sizes recently proposed for room-temperature 
superconductivity [To| . It is therefore unclear if a Meiss- 
ner effect could be observed via magnetometry experi- 
ments in single clusters, or below which cluster size it 
would become unobservable. What is clear, however, is 
that some clusters generally behave differently than their 
bulk counterpart. 

In cluster crystals or suspensions of clusters in an em- 
bedding medium, the situation is likely to be different 
from that of single clusters because of proximity effects. 
In Ref. [l9|, for example, the Meissner effect was ob- 
served for Pb nanoparticles embedded in an organic ma- 
trix. Alternative methods of detection which could be 
implemented by diamond magnetometry include mea- 
surements of the Knight shift or spin relaxation times, 
both of which are known to be indicators of the super- 
conducting phase transition. For example, Ti exhibits 
the Slichter-Hebel peak below T c [20]. Ti could also be 
used, in favorable cases, to report on the phonon spec- 
trum (see, e.g., Ref. [21]) and this would be useful to 
probe electron-phonon coupling mechanisms. 

The development of bulk room-temperature supercon- 
ductors out of atomic or molecular clusters may require 
a three-dimensional assembly of clusters into a lattice 
structure. Such three-dimensional assembly into cluster 
crystals has been demonstrated [22|. The study of indi- 
vidual clusters may be best performed by immobilizing 
individual clusters on a surface. Recently, "soft landing" 
of Ag56i clusters on a Ceo monolayer has been demon- 
strated [23j| . Another group has succeeded in immobi- 
lizing clusters on a surface Q by accurate placement 
and dimensional control in the assembly of clusters. This 
method would enable precise control of the number and 
position of clusters assembled on a surface, leading to 
a more systematic way to setup experiments that probe 
cluster-cluster interactions. 

Magnetic moments of atomic clusters are usually mea- 



sured in beam experiments using Stern-Gerlach deflec- 
tion [25]. An interesting variant would use surface mag- 
netic measurements to probe the cluster while it interacts 
with the surface or in the presence of an embedding ma- 
terial. Probing single clusters rather than populations 
would enable us to study heterogeneity in the proper- 
ties among an ensemble of clusters rather than measuring 
ensemble- averaged properties. In addition to supercon- 
ductivity and magnetism, one could also probe insulator- 
to-metal transitions as a function of cluster size and tem- 
perature. As atoms are added to the cluster, the number 
of delocalized electrons increases and the cluster's ini- 
tially paramagnetic state fluctuates and eventually be- 
comes diamagnetic. 

The magnetic properties (magnetization, M) of mag- 
netic materials are often measured with a susceptometer 
under an applied field (B). The M — B curve is typically 
probed using a SQUID (superconducting quantum inter- 
ference device) gradiometer oscillated sinusoidally near 
the sample and cycling the applied field. For clusters, a 
similar kind of experiment may be possible in which a 
bias field generates a strong enough magnetic polariza- 
tion M to observe the diamagnetism/paramagnetism of 
the cluster at the sensor's location. Let us consider the 
example of a Nag2 cluster (N = 92) placed in a magnetic 
field, B = 1 T. Assuming spherical clusters, the mag- 
netic moment of a diamagnetic cluster is given by the 
formula |2ffl 
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where R = ao^A/" 1 / 3 w 9.32 x 10" 10 m for Na 92 . r s is 
the Wigner-Seitz radius (r s = 3.9 for sodium) and ao is 
the Bohr radius, m and e are the mass and charge of the 
electron, respectively. This leads to ji^ia = —3.8 x 10 -25 
A-m 2 . If the sensor is placed at a distance of, for example, 
10 nm from the cluster, the induced field to be measured 
for the diamagnetic state is 
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or Bdia ~ 3.8 x 10 -5 T at 1 nm distance. 

In the paramagnetic state, we can apply the formula 
of Pauli paramagnetism for a free electron gas which 
gives [26| 
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where n in this formula is the number of atoms per unit 
volume, namely n — 0.97 mol/22.989 cm 3 for sodium 
meta, \±b is the Bohr magneton and the Fermi energy for 
Na is E F = k B T F = 3.24 eV. 

This gives a magnetization of M = XparaB ~ 0.7 x 
10 -5 T in a 1 T applied field. The corresponding dipole 
moment is 
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M ■ V « 1.94 x 10 _26 A • m 2 



(4) 
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(radius of Nag2 cluster is 9 x 10 10 m) . The corresponding 
dipole field is 
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at 10 nm separation distance, or B w 2x 10 -6 T at 
1 nm distance. To measure a field of magnitude 10 _9 T 
with a sensor whose noise floor is 4 nT /y/Hz requires 
about 8 seconds of averaging. We note that although the 
paramagnetism is less than the diamagnetism, the phase 
transition is predicted to be accompanied by a shift from 
paramagnetic to diamagnetic states upon crossing T c . 
This is only true of clusters. Indeed, diamagnetism in 
metal clusters is different from bulk metals. A Meiss- 
ner effect may be observable in clusters with incomplete 
shells, due to paramagnetism at high temperatures [9j, 
and the transition to the superconducting state should 
be accompanied by the paramagnetic-diamagnetic tran- 
sition in the cluster moment. 

A number of practical issues arise for NV-based sensors 
when a large (1 T) external field is applied. The three 
main considerations are: 1) the high microwave frequen- 
cies required, 2) the field drift and 3) the background 
magnetism of the substrate. The first element has been 
successfully addressed by Neumann et al. [27j who have 
interrogated NV centers in diamond under a bias field of 
0.65 T using ~10 GHz microwave fields. The field drift 
can be mitigated by avoiding fields that are too high and 
unstable. The 0.65 T field of [27| is sufficiently stable 
for optically-detected magnetic resonance (ODMR) mea- 
surements. The background magnetism of the substrate 
can be mitigated by translating the measuring tip near 
the cluster and performing subtraction of the signal from 
neighboring points (see Fig. QJ3) [14]. This technique 
has been successfully implemented by Moler's group in 
a scanning probe micro-SQUID susceptometer [28]. Fi- 
nally, we note that the energy level structure of the NV 
center is fundamentally different at high fields (> 0.1 T). 
However, the level structure itself should not present a 
problem for magnetometry. 

A possible experimental setup for interrogating atomic 
clusters and cluster crystals is shown in Fig. A 
nanoscale diamond tip is positioned near the surface of 
the cluster crystal or immobilized clusters and the tip is 
scanned along the surface to provide spatially-resolved 
information. The magnetic field gradient can be mea- 
sured by scanning the tip across the surface and perform- 
ing background subtractions using neighboring points 
(Fig. HP). An equivalent scheme is to translate the sam- 
ple rather than than tip. Nanodiamonds with a single NV 
center could be used, but so could nanocrystals contain- 
ing multiple NV centers. For detecting superconducting 
phase transitions, we envisage an experiment in which 
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FIG. 1. (A) Experimental setup for surface diamond magne- 
tometry of clusters. The diamond is scanned across the sur- 
face and surface magnetism is measured by subtracting the 
background signal from neighboring points (B). (Or equiva- 
lent ly, the sample can be translated rather than the tip.) 



both the external field and the sample temperature are 
cycled. 

As the cluster size increases the properties approach 
those of the bulk. A bulk material consisting of embed- 
ded superconducting clusters may exhibit properties that 
are closer to those of a bulk superconducting material in 
the sense that diamagnetic effects can be enhanced com- 
pared to those of the individual clusters. Strong fields 
may not be needed to generate an observable magnetic 
moment if diamagnetism is created by persistent cur- 
rents. In this case, the Meissner effect can be detected 
using a small external field. 

We now look at ways to measure the Meissner diamag- 
netism in a weak field using a diamond magnetometer. 
The experimental setup is shown in Figure El Green light 
(532 nm) from a 200-mW continuous- wave solid state 
laser is incident normal to the surface of a [lll]-oriented 
diamond crystal of dimensions 2 mm x 2 mm x 1 mm. 
The diamond was fabricated commercially by Sumitomo 
by the high-pressure, high-temperature (HPHT) method 
with an initial concentration of nitrogen impurities of less 
than 100 ppm, irradiated with a dose of 10 19 cm -2 , 3.0- 
MeV electrons, then annealed for 2 hours at 700 °C [29j . 
The fluorescence emitted by the diamond is recorded us- 
ing a Si photodiode. A 675 ± 75 nm interference fil- 
ter is used to attenuate the green light backscattered 
towards the photodiode. The BSCCO sample (BSCCO- 
2223, Bi2Sr2Ca2Cu30io), a 2.8 cm-diameter, 4 mm thick 
disc, is placed at a distance of appoximately 5 mm behind 
the diamond and parallel to the 2 mm x 2 mm face of the 
diamond. This particular sample-diamond geometric ar- 
rangement results in the largest magnetic-field shift due 
to the Meissner effect. An ambient bias field was gen- 
erated by a coil that surrounded the experimental setup 
and produced ~14 G at the diamond's location. 

The diamond was glued to a thin glass slide which 
was mounted on a printed circuit board with microwave 
striplines. A 200-micron diameter copper wire ran in 
between striplines ~1 mm from the surface of the dia- 
mond and delivered cw microwave field (~13 dBm power, 
2.895 to 2.915 GHz frequency- sweep range, 0.20-s sweep 
repetition time) from a Hewlett Packard 8350B Sweep 
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FIG. 2. Schematic diagram of the experimental setup. 



Oscillator with a HP 83525B RF Plug-In module. The 
photocurrent from the Si photodiode was recorded using 
an oscilloscope triggered upon each microwave frequency 
sweep. A Labview interface was used to download the 
data from the oscilloscope and record thermocouple volt- 
ages at regular intervals. Two Type-T thermocouples 
were used to measure temperature. The first thermo- 
couple was placed within 1 mm of the diamond. The 
second was held in close contact with the BSCCO sam- 
ple, and silicon-based heatsink compound was used to 
improve thermal contact between the thermocouple and 
the superconductor. 

Plots of the photocurrent vs. microwave frequency are 
shown in Fig.[3]for various BSCCO sample temperatures. 
For each such graph, the \ms = 0) O \ms = 1) transi- 
tion frequency from the [lll]-oriented NV centers was 
extracted by fitting a set of three Lorentzian lines sepa- 
rated by 2.2 MHz corresponding to the hyperfine splitting 
due to the 14 N nuclei (spin 7 = 1), which can be found 
in three possible orientations \mj = 0), \mj = ±1) with 
respect to the NV axis [29]. 

The \ms = 0) O \ms = 1) transition center frequency 
provides a measure of the projection of the local magnetic 
field B along the N-V axis, which causes a frequency 
shift gfisS • where 5 = 1 (spin), the magnetogyric 
ratio g\iB is ~2.8 MHz/G for the NV center electronic 
triplet spin state. This Zeeman splitting is in addition to 
the zero-field splitting (ZFS) H cf = D[S 2 Z - (1/3)5(5 + 
1)] with D w 2.87 GHz along the axis of the NV center 
(spin 5 = 1). For a superconductor below T c placed in 
a bias field, we expect the Meissner effect to alter the 
magnetic field around it, due to the exclusion effect of 
the superconductor. 

Our experiment consists of cooling the BSCCO sam- 
ple in liquid nitrogen for 30 s, then rapidly positioning it 
in close proximity (5 mm) to the diamond crystal, with 
the disc parallel to the crystal plane. The disc's axis is 
parallel to B. As the BSCCO sample warms up to room- 
temperature, measurements of the \ms = 0) <o> \ms = 1) 
transition in diamond is measured at regular intervals 
(every 4 s). With BSCCO-2223, T c is 105 K [30], and 



FIG. 3. NV-diamond detected microwave spectroscopy of the 
\ms = 0) <B> \ms = 1) transition for different BSCCO tem- 
peratures. Three peaks can be resolved which correspond to 
the hyperfine splitting from the 14 N nucleus. Depending on 
the temperature of the BSCCO superconductor, the center of 
the \ms — 0) \ms — 1) transition changes due to two com- 
peting effects: the Zeeman shift associated with the Meissner 
effect and a frequency shift associated with the diamond tem- 
perature, which may change as a result of the proximity to 
the BSCCO. See text for details. 



this gives us nearly 20 s at ambient temperature to ac- 
quire measurements in the presence of the Meissner effect 
before the phase transition is crossed (T > T c ) and the 
superconductive state is lost. The results of four differ- 
ent runs performed under similar conditions are shown 
in Figured! This protocol worked well for BSCCO, but 
for other materials, a better control of temperature, e.g. 
using a cryostat, may be required. 

The inset of Figure [H shows the least squares fit of the 
data to a sigmoid function 
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for one particular run. The free parameters of the fit are 
T c , /o and f\. T c indicates the phase transition, fo and 
/i are vertical offset and scaling parameters, respectively. 
The results for four runs are shown in Table H From the 
fitted value of T c , we determine a standard deviation of 
5.7 K. The average value of T c is 102 K. 

Run h [MHz] T c [K] r [K] f [MHz] 



1.3746 
1.5251 
1.0883 
1.4003 



101.3 
102.1 
104.3 
91.4 



4.53 2905.3 

4.20 2905.2 

4.71 2905.7 

5.33 2905.2 



TABLE I. Results of the least squares fit of a sigmoid to find 
the critical temperature. From these data, we find T c of (102 
± 3) K. 

Control experiments were performed with a Teflon 
disc cut to the same dimensions as the BSCCO sam- 
ple and cooled to liquid nitrogen temperature. Results 
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FIG. 4. BSCCO-2223 shows a phase transition to the su- 
perconducting state at T c of (102 ± 3) K. The four symbol 
types represent repeated runs on the same sample perfomed 
at different times. 



(not shown) did not exhibit any signs of a phase transi- 
tion up to room temperature, as one would expect of a 
non-superconducting material. 

The data in Figure 3] show a typical spread in the 
\ms = 0) <-> \ms = 1) transition frequency measure- 
ments found in our experiments. This spread is seen to 
be on the order of ca. 50-100 kHz in the high temper- 
ature range. Part of this drift is caused by drifts of the 
diamond temperature toward colder temperatures, which 
we estimate to be at most ca. 10 K, as the cold sample 
is in close proximity of the diamond. 

While performing the experiments in this study, we 
have noticed a statistically significant dependence of the 
ZFS parameter D of the NV center with temperature. 
This indicates that extra care should be taken in vari- 
able temperature studies, for the effects of dD/dT are 
otherwise indistinguishable from changes in the exter- 
nal field in the scheme we have used here. The problem 
could be mitigated or eliminated either by subtraction 
of the effects of a temperature drift from independent 
measurements of temperature, or by using gradiometer 
schemes. We have subsequently investigated the tem- 
perature dependence in details [3l| and found a value 



of dD/dT = -74.2(7) kHz/K in the temperature range 
280-330 K. In particular, we proposed another method to 
eliminate the effects of temperature drift of the ZFS co- 
efficients based on pulsed ESR (electron spin resonance). 

In conclusion, we have demonstrated the use of NV- 
diamond magnetometry for probing Meissner effects in 
a bulk superconductor and identified some issues associ- 
ated with such experiments. When positioned in close 
proximity to the surface of a superconductor, changes in 
the surface field can be detected with temperature. Care 
must be taken to ensure that the temperature depen- 
dence of the ZFS parameters does not adversely affect 
the measurement. 

We also estimate that it should be possible to probe 
magnetism in individual metal clusters. Nanoscale dia- 
monds with NV color centers present an opportunity for 
the study of magnetism and superconductivity in nano- 
materials, where sufficient resolution and sensitivity to 
observe properties of individual clusters could be ob- 
tained. In bulk superconductors, nanoscale NV-diamond 
magnetometry could be used to image vortices, param- 
agnetism and nucleation effects and could present advan- 
tages over micro- SQUID scanning techniques [32[ by en- 
abling variable temperature and hysteresis studies while 
maintaining close proximity to the sample throughout 
the cycle. Most importantly, high temperatures are in- 
accessible to SQUIDs without the use of thermal insula- 
tion layers that would prevent short-distance operation. 
While nucleation effects in the paramagnetic Meissner ef- 
fect are known to depend on cooling rate [32] and most 
likely exhibit hysteresis effects, magnetic imaging studies 
have not been carried out as a function of temperature 
across the phase transition. Such studies could shed light 
on the origins of the paramagnetism. Finally, the study 
of individual metal clusters and their interaction with en- 
vironment and surfaces could shed light on their roles in 
chemical catalysis. 
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